home *** CD-ROM | disk | FTP | other *** search
/ CU Amiga Super CD-ROM 19 / CU Amiga Magazine's Super CD-ROM 19 (1998)(EMAP Images)(GB)[!][issue 1998-02].iso / CUCD / Programming / LEDA / source / src / plane / _point.c < prev    next >
Encoding:
C/C++ Source or Header  |  1994-11-16  |  2.7 KB  |  132 lines
  1. /*******************************************************************************
  2. +
  3. +  LEDA  3.1c
  4. +
  5. +
  6. +  _point.c
  7. +
  8. +
  9. +  Copyright (c) 1994  by  Max-Planck-Institut fuer Informatik
  10. +  Im Stadtwald, 6600 Saarbruecken, FRG     
  11. +  All rights reserved.
  12. + 
  13. *******************************************************************************/
  14.  
  15.  
  16.  
  17. #include <LEDA/segment.h>
  18. #include <math.h>
  19. #include <ctype.h>
  20.  
  21. static const double eps = 1e-10;
  22.  
  23.  
  24. //------------------------------------------------------------------------------
  25. // points 
  26. //------------------------------------------------------------------------------
  27.  
  28.  
  29. point_rep::point_rep()  { count=1; x = y = 0.0; }
  30.  
  31. point_rep::point_rep(double a, double b) 
  32. { x = a; 
  33.   y = b; 
  34.   count = 1; 
  35. }
  36.  
  37.  
  38.  
  39. point::point()                  { PTR = new point_rep; }
  40. point::point(double x, double y){ PTR = new point_rep(x,y); }
  41. point::point(vector v)          { PTR = new point_rep(v[0], v[1]); }
  42.  
  43. double point::angle(const point& q, const point& r) const
  44. {
  45.   double cosfi,fi,norm;
  46.   
  47.   double dx  = q.ptr()->x - ptr()->x; 
  48.   double dy  = q.ptr()->y - ptr()->y; 
  49.  
  50.   double dxs = r.ptr()->x - q.ptr()->x; 
  51.   double dys = r.ptr()->y - q.ptr()->y; 
  52.   
  53.   cosfi=dx*dxs+dy*dys;
  54.   
  55.   norm=(dx*dx+dy*dy)*(dxs*dxs+dys*dys);
  56.  
  57.   cosfi /= sqrt( norm );
  58.  
  59.   if (cosfi >=  1.0 ) return 0;
  60.   if (cosfi <= -1.0 ) return M_PI;
  61.   
  62.   fi=acos(cosfi);
  63.  
  64.   if (dx*dys-dy*dxs>0) return fi;
  65.  
  66.   return -fi;
  67. }
  68.   
  69.  
  70. point point::rotate(const point& origin, double alpha) const
  71. { if (origin == *this) return *this;
  72.   segment s(origin,*this);
  73.   return s.rotate(alpha).end();
  74. }
  75.  
  76. point point::rotate(double alpha) const
  77. { return rotate(point(0,0),alpha); }
  78.  
  79. point point::translate(double alpha, double d) const
  80. { double dx = cos(alpha) * d;
  81.   double dy = sin(alpha) * d;
  82.   return point(ptr()->x+dx,ptr()->y+dy);
  83.  }
  84.  
  85. point point::translate(const vector& v) const
  86. { return point(ptr()->x+v[0],ptr()->y+v[1]);
  87.  }
  88.  
  89. double point::distance(const point& p)  const
  90. { return hypot(p.ptr()->x - ptr()->x, p.ptr()->y - ptr()->y); } 
  91.  
  92. double point::distance() const
  93. { return distance(point(0,0)); }
  94.  
  95.  
  96. int point::operator==(const point& p) const 
  97. { return (fabs(ptr()->x - p.ptr()->x) < eps && fabs(ptr()->y - p.ptr()->y) < eps); }
  98.  
  99.  
  100. ostream& operator<<(ostream& out, const point& p)
  101. { out << "(" << p.xcoord() << "," << p.ycoord() << ")";
  102.   return out;
  103.  } 
  104.  
  105. istream& operator>>(istream& in, point& p) 
  106. { // syntax: {(} x {,} y {)}
  107.  
  108.   double x,y; 
  109.   char c;
  110.  
  111.   do in.get(c); while (in && isspace(c));
  112.  
  113.   if (!in) return in;
  114.  
  115.   if (c != '(') in.putback(c);
  116.  
  117.   in >> x;
  118.  
  119.   do in.get(c); while (isspace(c));
  120.   if (c != ',') in.putback(c);
  121.  
  122.   in >> y; 
  123.  
  124.   do in.get(c); while (c == ' ');
  125.   if (c != ')') in.putback(c);
  126.  
  127.   p = point(x,y); 
  128.   return in; 
  129.  
  130.  } 
  131.  
  132.